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Abstract 



O ] We study the self-interaction effects for the Dirac particle moving in an external field created 

Q by static charges in (1 + 1)- dimensions. Assuming that the total electric charge of the system 

"^ vanishes, we show that the asymptotically linearly rising part of the external potential responsible 

+-' for nonexistence of bound states in the external field problem without self-interaction is cancelled by 

P-ithe self-potential of the zero mode of the Dirac particle charge density. We derive the Dirac equation 

r-| which includes the self-potential of the non-zero modes and is nonlinear. We solve the spectrum 

IJ problem in the case of two external positive charges of the same value and prove that the Dirac 

k> particle and external charges are confined in a stable system. 
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1 Introduction 

The problem of motion of the Dirac particle in an external field can be studied in different approx- 
imations. In the vanishing self-interaction approximation we neglect the effects of self-interaction, 
i.e. of the fact that the Dirac particle creates its own radiation field and interacts with it. Only the 
interaction between the Dirac particle and external field is taken into account. If the external field 
is created by static charges, then this interaction is described by the Coulomb potential. 

However, the vanishing self-interaction approximation cannot be applied universally; in some cases 
it leads to incorrect results and even to paradoxes. In particular, if we take in (1 + 1)- dimensions a 
system of two spin-1/2 particles of opposite charges coupled by instantaneous Coulomb interaction 
and assume that one of the particles is much heavier (proton), while the lighter particle (electron) 
moves in the Coulomb field of the heavier one, then it turns out that the system has not discrete 
energy levels [1]. In other words, in (1 + l)-dimensions the hydrogen atom does not exist. That 
happens not only for hydrogenlike systems with an infinitely heavy source of potential, but also for 
positroniumlike systems. 

It is therefore of principal importance to include self-interaction. We need the self-field effects to 
obtain the full picture of the interaction between the electromagnetic field and Dirac matter as well 
as to make that picture selfconsistent. In the self- field formulation, the electromagnetic field has as 
its source all the charged particles which in turn move in this field. The total electromagnetic field 
is a sum of an external and a self-field parts. The external part is created by some external sources 
which are not dynamically relevant, and the self-field is created by the Dirac particle itself. 

In the present paper we aim to study the effects of self-interaction for the Dirac particle in (1 + 1)- 
dimensions in connection with the result of [1]. We want to determine whether in the presence of the 
self-field the Dirac particle and external charges can be confined in a stable system characterized by 
discrete energy levels. 

Models in (1 + l)-dimensions are interesting as simpler models for discussion of different aspects 
of realistic particle physics models in (3 + l)-dimensions. At the same time, (1 + l)-dimensional 
models are interesting in their own right and have some pecularities which make physics in (1 + 1)- 
dimensions different in principle from one in (3 -|- l)-dimensions. One of these pecularities is that the 
Coulomb potential on line is linearly rising at spatial infinities. Just the linear Coulomb potential is 
responsible for the paradox mentioned above. 

One-dimensional models of spin-1/2 particles have applications in condensed matter physics, too. 
It is enough to mention the one-dimensional model of electronic liquid or the Thirring model. The 
problem of existence of bound states for a Dirac particle in the presence of a static charge distribution 
may be useful for understanding of formation and decay of bound states for these particles in other 
one-dimensional models with more complicated interaction. 

Our paper is organized as follows. In section 2, we neglect first self-interaction and consider the 
Dirac particle moving in a potential created by external static charges. We find asymptotics of the 
eigenfunctions of the Dirac Hamiltonian at spatial infinities. In accordance with [1], for any value 
of energy the eigenfunctions are not normalizable and cannot correspond to a discrete spectrum. In 
section 3, we introduce the self-field of the Dirac particle. Following the method of [2-4], we derive 
the Dirac equation which includes the nonlinear self-field term. In section 4, we study this equation 
for a specific choice of the external potential, namely for the case of two external positive charges 
of the same value. We solve the spectrum problem in the approximation when the discrete energy 
values are determined by the interaction between the Dirac particle and external field, while the 
self- interaction shifts these values by a small amount. Section 5 contains our conclusions. 



2 Dirac particle in external field 

The Dirac equation for a particle of charge e and mass m moving in an external field is 

[Yiificdf, - eA'^^') - mc^]^{x) = 0, (2.1) 

where {jj,, z/ = 0, 1), 7^" are 2x2 Dirac matrices, 

and the field iIj{x) = ilj{xi,t) is two-component Dirac spinor. The partial derivatives are defined as 
do = d/cdt, di = d/dxi. 

We assume that the external field is defined by 

A^^*(xi,t) = (Ag^* = A',^\x^),Ar' = 0), 

which is a time-independent scalar potential. We assume next that the potential is created by static 
charges, so we can use the following ansatz 

A-*(xi) = -^Q'^'\x,\ + A-*(a;i), (2.2) 

where Q^^^ is total external charge, and Aq^* is a potential which does not increase at spatial infinities, 
i.e. Aq^^{xi = ±00) are finite constants. 

In particular, one external charge q at the point xi = creates the current j'q^* = qS{xi) that 
corresponds to the scalar potential A^^^ = —■^q\xi\. Two external charges of the same value q taken 
at the points xi = a and Xi = —a create the current Jq^* = q{S{xi — a) + S{xi + a)) and potential of 
the form (2.2) with Q^^^ = 2q and 

g(|a;i| — a) for \xi\ < a, 
A^o"' = { (2.3) 

for |a;i| > a. 

We can easily check that the ansatz (2.2) is valid for an arbitrary number of external static charges 
gi, q2, qs, ... . Therefore, the potential for these charges can be always given as the sum of two parts, 
linearly rising and finite at spatial infinities. 

If we act on both sides of (2.1) by [—'-/^{ihcdfj, — eA^^^) — mc^], then we come to the second-order 
differential equation for ip [5,6] 



Tn c 



^ = 0, (2.4) 



where 

This equation looks like the Klein-Gordon one, except the additional term eS^^F^f^ with 



s^^^ ^Uyi\^^]_, 



In (3 + l)-dimensions, the spatial components of S^'^ are related to spin of the Dirac particle, so the 
additional term describes interaction between the particle spin and the electromagnetic field. 

In (1 + l)-dimensions, S*^^ has no spatial components, and we can not therefore introduce spin. 
The only nonvanishing component S^^ = ia, where a = 7^ = '~f^'~f^, allows us to introduce chirality. 
If we write ip in the component form as 

and define the operators 



then 



^^V'-U- ■ ^-^'-u. 



of Ff !^* is 



\ext 

T-iext '^^^0 c^^i 

-Tni — 7^ = O 



i.e. the upper component is of positive chirality and the lower one - of negative chirality. 

Moreover, there is no magnetic field in (1 + l)-dimensions, so the only nonvanishing component 

Text 

where S'^^^ is the external electric field. The additional term in (4) becomes 

indicating that positive and negative chirality components are coupled differently to the external 
electric field. 

In the Hamiltonian form the equation (2.1) reads 

d f d \ 

Hij = ih—iJ = ihca—- + (3mc^ + eA"^* ^, (2.5) 

at \ oxi J 

where (3 = 7°. The eigenvalue problem for the Dirac Hamiltonian, Hip = Eip, reduces to the problem 
of solving the system of two equations 

in— 1 Ml = —mcu2, (2.6a) 

( ^ d eA'^r?'^ -E\ ,^^, , 

—iTij^ 1 I U2 = —mcui. (2.Db) 

We easily decouple ui and U2 and rewrite these equations equivalently as 

Ml(2) =m^C^Ml(2), (2.7) 



s^^±A...+ ^-4?*-^V 



dxl c ' \ c 

the sign (+) corresponding to the positive chirality component, while (— ) to the negative one 



If the Dirac particle and external charges are confined in a stable system, then equations (2.7) 
must reveal a set of bound states. For a discrete set of energies in the band \E\ < mc^, these 
equations must have solutions which decrease exponentially at infinities and are normalizable. 

However, for Q^^^ ^ 0, the external potential is asymptotically linearly rising at spatial infinities, 
so the term (e^((5^^*)^a:^)/(4c^) dominates in the equations (2.7) and prevents any bound state. 
Indeed, the asymptotics of Mi(2) for all possible energies in the band are determined by the equation 



dxl 



ext\2 



+ 



Ac^r 



~ •X' -\ 



Ml(2)(kl| 



OO 



0. 



(2. 



This is the inverted oscillator equation [7]. Its general solution can be expressed in terms of parabolic 
cylinder functions. One way to choose two linearly independent solutions of (2.8) is to take the real 
functions W^(0,xi) and M^(0, — xi) (we follow the notations of [7]) whose asymptotic behaviour is 
well-known: 



1^(0, xi -^ +oo) 
1^(0, a;i -^ -oo) 



Xi 

1 



Xi\ 






We can represent 'Ui(2)(|a;i| -^ oo) as linear combinations of W{Q^xi -^ ±oo) with arbitrary 
coefficients. Regardless of the choice for these coefficients, the normalization integral 



x-i 



dx\{\ui{x\)\^ + \u2{x'^)\^) 

diverges for large Xi. This means that the probability of finding the Dirac particle and external 
charges infinitely separated remains finite at all times. Consequently, the equations (2.7) have solu- 
tions for an arbitrary energy, but these solutions can never represent a bound state. 

3 Self-interaction 



Let us now introduce the self-field of the Dirac particle and investigate how it infiuences the bound 
state spectrum. To derive an equation for the Dirac particle in the presence of both the external and 
its own radiation fields , we start first with the action 



S 



dt 



dx,[i>{xi,t){Ythcd^ - mc^)iP{xi,t) - j';\xi,t)A''{xi,t) 



-iF,^(a;i,t)F^^(a;i,t)], 



(3.1) 



■ext 



where 

is a total current which includes both the Dirac matter current j^ = eip'-f^ip and the one created by 
external static charges. 

If we expand the field ip into the Fourier integral 



^{xi,t) 



1 



/27r 



dp-e^P'''^{p,t), 



then the charge density of the Dirac matter becomes 



jo{xi,t) = eij*{xi, t)^(xi, t) = — / dk- e-n>'^'p{k, t), (3.2) 

where 

p{k,t) = dpij*{p + k,t)ij{p,t). 

The zero momentum component of the density determines the matter charge 

/oo 
dxJQ{xi,t) = ehp{0,t). 
-oo 

Separating in (3.2) p{0,t) and p{k,t) with the non-zero momentums A; 7^ 0, we rewrite (3.2) as 

joixi, t) = Q5{xi) + joixi, t), (3.3) 

the density 

JQ{xut) = — dk-e-n^'^'pik^t), 

ZTT J-00 

p(A;,t)=p(A;,t)-p(0,t), 

corresponding to a zero charge, /^ dxij^^Xi, t) = 0. 

With the choice of the gauge (9^A^ = , the Maxwell equations take the form 

OA,=f;\ (3.4) 

and are solved by 

/oo roo 

dx[ / dt'D\xi -x[;t- t')f°\x\,t') 
-00 J —00 

= <^(xi,t) + A-*(xi,t), (3.5) 

where D'^{xi,t) is the causal Green's function 

aD%xi,t) = 6{xi)S{ct), 

D%x^,t) = -— — / dqo / dq,- ^—e-^^^otg^.i^i (g.g) 

(27r)^ J-00 J-00 ^o ~ Q'l + ^f 

and □ = d^d'^. 

Substituting (3.3) and (3.6) into the expression for the electromagnetic field, we get its external 
and self-field parts in the form 

A'';\x^,t) = (^-iQ-*|xi| + A-*(xi)) 6,0, 



that agrees with (2.2), and 



Aff{x,,t) = -^Q\x,\S,o + Aff{xut), 



with _ _ 



/oo roo _ 

dx[ / dt'D^xi -x[;t- t')j^{x[,t'), 
-00 J —00 



and 

Putting both parts together, we see that if the total charge of the system vanishes 

Qtot ^ g + Qext ^ Q^ (3 7) 

then the hnearly rising potentials produced by the Dirac particle and external charges cancel each 
other, and for the total electromagnetic field we get 

A^ix^,t) = Aff{x^,t) + A',^\x,)5^o- (3.8) 

One of the Maxwell equations (3.4) is the Gauss' law 

^f^.r. (3.9) 

where the electric field also consists of two parts, self-field and external, 

^ = Foi = S'^^^ + S^^'K 
Integrating the equation (3.9) over Xi, we obtain 



i.e. the difference in the values of the self-field electric field at the ends of Xi-line balances the one 
of the external electric field, and for the total electric field S{+oo) = S{—oo). 

The vanishing of the total charge allows therefore to balance the sources of electric flux and is 
important from a physical point of view. For Q'^^^ ^ 0, the balance is destroyed, and, as we have 
seen in the previous section, the Dirac particle and external charges escape from each other and can 
not create a stable system. 

The condition (3.7) has its analogues in the second- quantized version of different (l+l)-dimensional 
models. In the Schwinger model [8] , the total electric charge is known to be zero on the physical states 
[9,10]. Because of Schwinger [8], when a charge is inserted into the vacuum, the accompanying electric 
field polarizes the vacuum producing complete compensation of the charge. In (1 + l)-dimensional 
QCD, we restrict ourselves to color neutral states, since the presence of uncompensated color charge 
in space leads to a growth of the fields at infinities and makes the total energy of the system infinite 
[11]. 

Inserting (3.8) into the action and using a partial integration, we have 

S = So + Sself, 

/OO /.QO _ _ _ 

dt / dxi['ilj{xi,t){Yi^cdf, - mc^)tlj{xi,t) - etlj{xi,t)-f^'tlj{xi,t) ■ A""^^ 
-OO J —OO 



\fe.tATi (3.10) 



dt / dx,[-ei^{x,, t)rH^,,t)Aff - -Flffp^^^f 

-OO J —OO 4 



with F^^ = dxAf^^^ — d^A'l^ and A'jf* = Aq^^S^o- Variation of this action with respect to the Dirac 
field yields the following Dirac equation 



^{■i^cdf, - eA"^^') - mc^\ ^{x) = e^Af-fipix). (3.11) 



It is essential that neither the hnearly rising part of the external electric field nor the one of the 
self-field enter this equation. 

Since AJf''-^ is expressed in terms of the current j^, (3.11) is a non-linear integral equation for 
ip. In the next section, we will continue our study of the equation (3.11) for a specific choice of the 
external electric field. 

4 Example 

1. As an example, let us consider in detail the case of two external charges of the same value g > 
at the points Xi = ±a and find first the spectrum for the vanishing ^ff'-'^. 

With the external potential Aq^* given by (2.3) (see Fig. 1.), the Dirac Hamiltonian is invariant 
under a modified parity transformation generated by 

f ^ pi, (4.1) 

where the action of / on any function of Xi is defined as 

ifixi) = fi-xi). 

The Hamiltonian eigenfunctions must be therefore even or odd under parity reversal. Since the 
transformation for ip is 

ij{xi, t) -^ tp'{xi,t) = P^/j{-Xi,t), 

we get 

Ui{Xi,t) =U2{-Xi,t) 

for even eigenfunctions, and 

Ui{Xi,t) = -U2{-Xi,t) 

for odd ones. 

To solve the system (2.6a-b) with A^^^, we have to find general solutions in the regions —a < 
Xi < a, Xi > a and xi < —a, and match them at the boundaries so that the eigenfunction and its 
first derivative is continuous. We want also to distinguish the regions of positive and negative values 
of Xi, and so introduce 

^1(2) = wi(2)(a;i > 0), 
^1(2) - ^1(2) (a^i < 0). 

To simplify the notation, we define in the region < xi < a the new dimensionless variable 

z^ = e-'^*.[^{xi-a-—). (4.2) 

V he eq 



The system (2.6a-b) then becomes 



' d \ 2 ^ 

■^ + z+\ui = iJ ^e*^M+, (4.3a) 



d \ / 2 TT 

+ z+]u^ = i\ —e'^uf, (4.3b) 



dz^ ' ^ ' V A 



where A = {2eqh)/{m'^c^) is a diinensionless constant. Decoupling uf and U2, we obtain the system 
of two second-order differential equations 



92 



Uw)\ + 



5^2 "1(2) 



■'+ 



{2v±l)- zl 



«1(2) = 



with V = —i/A. 

The system (4.4) is solved by hypergeometric functions. If we introduce 



(4.4) 



'1(2) 



- ^+5 

6 2+ ■ U 



1(2)' 



then (4.4) reduces to 



d^ 



l/.^^^i,. + (2-y. 



y+ 



dy- 



d 
dy+ 

d 



^1, 



— Vi 

2 ' 



0, 



-2,. + (2-y+)^-2,.+ 



Z/-1 



'-V2,+ = 0, 



(4.5a) 



(4.5b) 



which are just the hypergeometric function equations. For the first equation, the linearly independent 



solutions are F{—^; i; y^) and yj F^^-^-; |; |/_|_), while for the second one F{K^] i; |/+) and y/ F(l 



1/2) 



2! 2^?/+^ 



A similar analysis can be performed in the region —a < a;i < 0. With the variable 

eq' 



z_=e 4./-— (xi + aH ), 

\ nc 



the equations for -^^.2) I'sad 
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The linearly independent solutions for v^ 



(2z/^l) -zi 



^1(2) = 0- 



1/2 r./!-,.. 3. 



(2) 



e^^- ■ iz-(2) are F(if^; i;|/_), ^^^(1 - 



v_. 3. 

2' 2' 



(4.6) 
yj) and 



-^(~f 5 2' ^-)' ^- ^i^^'i 2' ^-)' correspondingly, where |/_ = -2^. 

Taking linear combinations of these solutions, we can construct eigenfunctions 

,± \ 



r^ 



Ml 



which fulfil the matching conditions 



i}^{xi = +0) = ij-{xi = -0), 



dxi 



{xi = +0) 



dtp 
dxi 



{xi = -0), 



and are even or odd under parity reversal. 

The even eigenfunctions up to a constant factor are 



i^i 



e 2^+ 



^e^fF(-|;i;4) + G^+F( 



l-u . 3. 2 A 

2 '2' +/ 



z^F{l - f ; I; 4) - G Jf e^f ^(if^; 1; 4) 



(4.7a) 



for < xi < a, and 



1,2 



-z_F{l - i^; i; zl) - Gj^e^iFiK^; i; z^ 






2z/^,F(l - I; f ; 2//^^) + F(-|; i; 2//^^ 



for —a < a;i < 0. The constant 

2,yz F(l — 

G(A,..;_ 2^,„F(if^;|;2z..^)_F(i^;i;2z..^) 

is modulo 1, \G\^ = 1, and Za = {E + eqa)/{mc^). 
The odd eigenfunctions are 



1,2 






,/,+ _ „-fz+ I V 2 "" -" V 2' 2' ^+/ ' ^^+-^ V 2 ' 2' ^+V I (. r. \ 



for < Xi < a, and 



^_F(1 - f ; I; zl) + G^f e^f F( V; ^; z' 
/|e^fF(_|;i;^2) + G^_^(i^.|.^2 



^-=^"^^M ~/i.,"^:;. 2/k .i-.:i' ^^ i (^.sb) 



for —a < a;i < 0, where 

For xi> a ii can be checked that 



with 



^i™=»r,-rU--- (4.9) 



K = \J'm?-c^ — E"^ 



satisfies the Dirac equation. Matching (4.7a) and (4.9) at Xi = a and ehminating s we obtain the 
equation that determines the spectrum of the even bound states 



E = mc^ cos {Xg{E)) (4.10) 

where 

Xg = arg 



2uGzoFC-^; |; 2uzI) - F(-|; |; 2uzl] 



y2uzoF{l - |; |; 2iyzi) + GFC-^; 1; 2uzi) J 
and zq = E/{mc^). The matching condition at xi = —a for (4.7b) and 

E_^^y~^"''\ xi<-a, (4.11) 

gives the same spectrum equation. 

In a similar way we can derive the equation that determines the spectrum of the odd bound 

states: 

E = mc^ cos {\g{E)). (4.12) 

10 



2. With the self-field term A'1^''-^ the Dirac equation is nonlinear in ip, so the spectrum problem 
becomes essentially more complicated. As in [2-4], we can solve the problem in the approximation 
when the self-interaction contribution to the energy spectrum is very small with respect to the 
contribution of the interaction between the Dirac particle and external field. 

Let us assume that the equations (4.10) and (4.12) have iV+ and iV_ solutions , correspondingly, 
i.e. for A^f-I" = there are N. even and N_ odd bound states. The total number of discrete states in 
the band |-E| < mc^ is then A^ = A^^_ -|- A^_. Let us denote the normalized bound state eigenfunctions 
by V"^^*, n = 1,N. Up to a normalization factor ip^^* coincide with (4.7a-b) for even and with (4.8a-b) 
for odd states in the region —a < Xi < a. 

In (3 -|- l)-dimensions, the self- field effects are of order of the fine structure constant and higher. 
This allows us to restrict our calculations to the first order of this constant. In our study, to make the 
self- field effects small we assume that |e| <^ g, so the self-interaction shifts the bound state energies 
by a small amount 

E„ = Ef * + AE^*^'^ (4.13) 

and does not change the number of states. In (4.13) E^^^ are discrete spectrum energies of the external 
field problem without Af^''-^ ■ The eigenfunctions of the discrete spectrum cannot be characterized 
now by a definite parity, because the self-interaction term in the Dirac Hamiltonian is not in general 
invariant under parity reversal. 

For the nonlinear Dirac equation the superposition principle does not hold. It is not possible 
to expand the exact solutions of this nonlinear equation as a superposition of the linear equation 
solutions ipf^^^ with arbitrary time-dependent coefficients, even though ipf^^ form a complete set. 
What can we use now is the Fourier expansion in the time coordinate [2-4] 

^(xi, t) = Y. V^n(xi)e-^^"* + / dEi,{x^, E)e-^nEt 

n=l -^-"^ 

+ / dEij{xi,E)e-i^' (4.14) 

in which the time behaviour is known, and of the form exp(— |i?t). The functions ipnixi), ip{xi, E) 
and the energies En are unknown. 

To derive an information about the spectrum it is simpler and more convenient to work with the 
action rather than with the Dirac equation. Since manipulations which we are going to do below are 
valid in both discrete and continuous spectrums, we can use instead of (4.14) the following compact 
expression 

ij{x^,t) =Y.Jn{xi)e-^nE"\ (4.15) 

where J2n means summation over discrete states and integration over continuous ones. 

If we insert the Fourier expansion into the action, then up to the terms depending only on the 
external field we obtain 



) 



where ujnm = En — Em, and 

p^Q—— /"oo poo poo poo 

Sself = 77-J2r, mrJ ^W ^^M ^^ / ^^1 bnm(a;i)^^(a;i - x[; t - t')jrsA^l) 

I "'™'''* J-oo J -oo J -oo J -oo 

-jL(xi)/^^(0;t-t')jV.,o(:r;)]e^(--*+--*'), 
11 



where 

After time integration for 5*0 we find 

5*0 = 27r^V / dxiilj^{xi){Em-~ H)lljjn{Xl)6{uJnm), (4.16) 

with H given by (2.5) . If the ipn{xi) were solutions of the external field problem with the vanishing 
self-potential AJ^^^ , i.e. ■?/'^^*(a;i), then this expression would be zero for En = E'^^. The '4''^^{xi) 
minimize 5*0 alone. However, now the entire action S , of which 5*0 is only one term, must be 
minimized as a whole. 

Time integrations in Sg^if can be performed using (3.6) and we can write the self-field part of the 
action entirely in terms of the Fourier components of the currents 

p^fj^-—— roo (Iq 

^.e./ = ^^E„^,/(^n- + ^-) / ^-^^^-^--■{TUqi)TrsA-Ql)-TL{0)T^sm, (4-17) 
Z ! ) 1 J—OO „uy^^, c/i ~\~ It 



TrLiQi)^ rfxij^„,(xi)e^^^-^ (4.18) 



where 



The ipn{xi) here are still exact solutions for which the action S will vanish identically. 

In the approximation of small self-interaction contribution, we can solve the spectrum problem 
iteratively. To lowest order of iteration we replace ipn{xi) with %Ij^^[xi) and then solve for En which 
have the form (4.13). The number of discrete states without and with the self- field term A^^^^ 
is assumed to be the same, and transitions between discrete and continuous states are neglected. 
Omitting the integration over continuous states and using the orthonormality of ip'^^-, we write 5*0 in 
the form 

N 

So = 27rh J2 ^Ef^ ■ K^nm)5nm- (4.19) 

?i,m=l 

In the self-field part of the action we separate the terms according to En = Em, Er = Eg and 
according to Eg = En, Er = Em, the only two ways of satisfying the overall 5- function for discrete 
spectrum energies. And since S* = to this order of iteration we can solve for ^En^^ ■ The action 
and the total energy of the system are related by a 5-function. Cancelling this (5-function as well as 
the sum over n to obtain the energy shift of a fixed level n, we get 



'-^E = y^ / dqiV^ ■ Gnn,mm{(ll) Z — Z^ / 1 i'^l i [{Gnm,mn{(ll) 

47r _-, J-oo ai Svr _-, J-oo \Q^ \ -uJr,m — loi I 

(4.20) 



a fic _L 

+ Gm.n,nm{qi)) + i^— ^ — R-G 

(m<n) 



^nm 



C C 



where 



Gnm,rs{qi) = ^nm'''(Q'l)^ri^/i(~Q'l) "~ ^nm' (O)^rT ' (0)) 

V stands for the principal value prescription, while Re means the real part. The subscript ext here 
indicates that the Fourier components (4.18) must be calculated by using ipn^^. 

The last term in (4.20) contributes if m < n. This shows that only the ground state n = 1 of 
our iV-level system is stable. The energy shift for the excited states n = 2, ..., A^ is complex, so the 
imaginary part of the shift can be identified with the spontaneous emission of the excited states due 
to self-interaction. 
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5 Conclusions 

1. We have studied the bound state problem for the Dirac particle moving in both external and 
its own radiation fields in (1 + l)-dimensions. We have shown that if the total electric charge of 
the system vanishes, then the asymptotically linearly rising part of the external potential which was 
responsible for nonexistence of bound states in the external field problem without self-interaction is 
cancelled by the self-potential of the zero mode of the Dirac particle charge density. The resulting 
Dirac equation includes only that part of the external potential which is finite at spatial infinities 
and also the self-potential created by the non-zero modes of the charge density. 

We have proved that this equation has a set of solutions which show that the Dirac particle 
and external charges are confined in a stable system. These solutions correspond to energy levels in 
the band \E\ < mc^. Only the lowest level is precise. The higher levels have a nonzero linewidth 
which manifests itself as spontaneous emission. This picture is characteristic for hydrogenlike atoms. 
According to the self-field approach, the hydrogen atom has no precisely defined sharp energy levels, 
other than the ground state [2-4]. The excited states cannot be stable due to radiation reaction. 

2. We have solved the bound state problem for the external potential created by two positive 
static charges of the same value. All other cases when we have two or more external charges can be 
considered analogously. 

The exceptional case is one external charge at the origin. The potential created by this charge 
consists only of the asymptotically linearly rising part, so the corresponding nonlinear Dirac equation 
does not include any external field. This means that the approximation when the self-interaction 
contribution to the spectrum is small with respect to the one of the external field cannot be applied. 
We need here to look for other ways to solve the nonlinear Dirac equation. This problem remains 
open. 

3. Vanishing of the total electric charge is the condition necessary for the existence of bound 
states of the Dirac particle and external static charges on line. If the total electric charge is not zero, 
then the Dirac equation does not reveal bound state solutions. Therefore, in the Dirac theory on line 
bound states can be only neutral. This is one of the pecularities of one-dimensional physics. 

We can consider charged states as well, but these states are not stable. The presence of an 
uncompensated electric charge leads to the linearly rising Coulomb potential in the Dirac equation 
and to a potential of the inverted oscillator type in the corresponding Schrodinger equation. Such 
potentials are known to allow only metastable states [12]. So for limited times the Dirac particles 
and external charges can be confined in a charged metastable state. Metastable states can decay into 
stable ones. If particles emitted in decay take away a charge equal to the charge of the metastable 
state, then the remaining particles can create a neutral bound state. 
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Figure Captions 

Figure 1 

The potential Ao{xi) in tlie case of two external charges of the same value g > at the points 

Xi = ±a. 
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FIG. 1. 
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